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O ; Abstract 

The popular demonstration involving a permanent magnet falling through a conducting pipe 
is treated as an axially symmetric boundary value problem. Specifically, Maxwell's equations 
are solved for an axially symmetric magnet moving coaxially inside an infinitely long, conducting 
cylindrical shell of arbitrary thickness at nonrelativistic speeds. Analytic solutions for the fields are 
developed and used to derive the resulting drag force acting on the magnet in integral form. This 
treatment represents a significant improvement over existing models which idealize the problem as 

5/5 ■ a point dipole moving slowly inside a pipe of negligible thickness. It also provides a rigorous study 

of eddy currents under a broad range of conditions, and can be used for precision magnetic braking 
applications. The case of a uniformly magnetized cylindrical magnet is considered in detail, and 

^ \ a comprehensive analytical and numerical study of the properties of the drag force is presented 

for this geometry. Various limiting cases of interest involving the shape and speed of the magnet 
and the full range of conductivity and magnetic behavior of the pipe material are investigated and 

r*| ■ corresponding asymptotic formulas are developed. 
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I. INTRODUCTION 

A. Background and Significance 

This work is based on the popular experiment which demonstrates magnetic damping 
by means of a permanent magnet falling through a conducting pipe [1, 2, 3, 4]. Such 
an arrangement has long been a favorite for demonstrating such topics as Faraday's law 
of induction, Lenz's law, eddy currents, inductive heating and magnetic damping. The 
underlying physical process is induction heating: the moving magnet creates a changing 
magnetic flux in its vicinity which induces circulating eddy currents within the pipe wall, 
thereby causing ohmic dissipation and generating a drag force on the magnet by virtue of 
energy conservation. The force itself is of course manifested as the action of the magnetic 
field generated by the eddies on the permanent magnet, with its retarding nature understood 
as a manifestation of Lenz's law. A particularly vivid picture of this mechanism emerges 
if one views the magnet as an assembly of circulating atomic currents moving through the 
pipe. Lenz's law then implies that the induced eddies in the pipe wall counter-circulate 
ahead of the moving magnet and co-circulate behind it. But this implies that the moving 
magnet is repelled in front and attracted in rear, hence acted upon by a retarding force. 

The main contribution to the literature on this subject is Saslow's paper [2] in which he 
treated the problem of eddy currents in thin conducting sheets quite generally, and as an 
example provided an approximate calculation of the drag force on a magnet falling inside a 
conducting pipe. MacLatchy et al. [3] used other methods to derive Saslow's result for the 
terminal speed and pointed out a sizable discrepancy between the predictions of this result 
and their measured values. These authors introduced a numerical modeling of the magnet 
which significantly improved the agreement. Following Saslow's calculation [2], analytical 
treatments of the magnet-pipe problem have considered the magnet as a point dipole moving 
at low speeds and the pipe as infinitely thin-walled and long. These assumptions imply that 
the only significant length parameter in the problem is the interior diameter of the pipe, 
and lead to a simple expression for the drag force. We will refer to this limit as the idealized 
model and derive it as a limiting case of our general solution in appendix B. The main 
sources of inaccuracy in this model are the point-dipole and thin-wall assumptions, with the 
low-speed approximation a potential source of inaccuracy as well. 

B. Objectives and Limitations 

Our objective in this work is to develop a rigorous formulation of the magnet-pipe system 
that avoids the approximations of the idealized model. To that end, we treat the case of 
an axially symmetric permanent magnet moving coaxially inside an infinitely long, conduct- 
ing cylindrical shell of arbitrary thickness. By an axially symmetric magnet we mean a 
permanent magnet whose magnetization vector has an axially (or azimuthally) symmetric 
magnitude and a uniform direction parallel to the symmetry axis. Since any practical real- 
ization of this model will likely involve magnet speeds far smaller than the speed of light, we 
will restrict our treatment to nonrelativistic speeds, v/c <C 1, where v is the magnet speed 
and c is the speed of light. On the other hand, we are including the possibility of the magnet 
being projected into the pipe thus allowing much higher speeds than can be attained by a 
magnet that starts to fall from rest under gravity. The restriction to nonrelativistic speeds 
implies that we are dealing with quasi-static sources and fields where displacement currents 



can be neglected [5]. This is so not only in the interior of the pipe where the inclusion of 
the displacement currents would amount to a minute correction of the order of (v/c) , but 
also within the pipe wall where conduction currents dominate displacement currents. To 
provide a basis for the latter assertion, we note that the basic time scale generated by the 
motion of the magnet is of the order of Ri/v, where Ri is the inner radius of the pipe. This 
time scale corresponds to a frequency of the order of v/R\, which implies that the ratio of 
displacement to conduction currents is of the order of e^v/aRi where a is the conductivity 
of the pipe wall [6]. For common metals and with 1 < Ri < 10 cm, this ratio is in the 
3(10~ 7 — 10~ 10 )f/c range. Since v/c <C 1 by assumption, the ratio in question is seen to be 
vanishingly small. 

The assumption of an infinitely long pipe is unavoidable if a manageable solution is 
desired. The error resulting from this assumption, on the other hand, is small if the magnet 
is not close to the pipe ends. To provide a basis for this assertion, one can use the idealized 
solution, appendix B, to estimate the dissipated power within the pipe segment that extends 
from its actual end to its idealized end (which is infinitely far). The ratio of this quantity 
to the total dissipated power is then an estimate of the leading finite-length correction to 
the drag force. This ratio is found to be of the order of (R\/D) , where D is the distance 
from the magnet to the near end of the pipe and R\ is the pipe's interior diameter. For 
Ri/D < 1/4, for example, the expected correction would be of order 10 -4 which is quite small 
as claimed. Similarly, the quasi-static approximation implies that the radiated power from 
the magnet-pipe assembly is negligible. To get an idea of the magnitude of such radiation, 
consider a permanent magnet with a magnetic dipole moment m moving longitudinally (i.e., 
parallel to m) through free space at nonrelativistic speeds. The radiated power from such 
a source can be calculated and is found to be fio'm 2 v 2 /30'KC 5 where v and v are the first 
and second time derivatives of the magnet's speed. This quantity should be compared to 
the ohmic dissipation in the pipe wall which for this purpose may be estimated using the 
idealized model. Using the solution to the idealized model given in appendix B, we estimate 
the ohmic dissipation rate to be 45/igm 2 w 2 crs/1024i?i 4 where s is the thickness of the pipe 
wall. For reasonable values of the parameters, the radiated power is totally dominated by 
the ohmic dissipation rate hence confirming the expectation that radiation is completely 
negligible in this problem. 

The case of a uniformly magnetized cylindrical magnet is considered in detail, and a com- 
prehensive analytical and numerical study of the properties of the drag force is presented 
for this case. Various limiting cases of interest are explored and appropriate asymptotic for- 
mulas are developed. The results reported here are supplemented with a computer program 
posted on the web [7] which can be used to compute the drag force for this case. 

The rest of this paper is organized as follows: section II A deals with the characterization 
of the source currents of a permanent cylindrical magnet in axial motion, and II B with 
the vector potential of such a magnet in the quasi-static limit. The electromagnetic fields 
of the magnet-pipe system are found in section III and the magnetic drag force is derived 
IV. Section V presents the results of numerical computations of the drag for the case of a 
uniformly magnetized cylinder as well as a detailed discussion of its dependence on magnet 
shape and speed, and also on the material properties of the pipe. Concluding remarks 
are presented in section VI. Several limiting cases of the drag force are considered in the 
appendices and corresponding asymptotic formulas are derived. 



II. SOURCES AND FIELDS OF A MAGNET MOVING IN FREE SPACE 

The first step in our analysis is the calculation of the electromagnetic fields of an axially 
symmetric permanent magnet moving along its axis of symmetry in free space as input for 
the magnet-pipe configuration. 



A. Source currents of a moving magnet 

By a permanent (or hard) magnet is meant a ferromagnetic material whose magnetiza- 
tion does not change when immersed in external fields, electromagnetic or gravitational. 
In practice, this requirement is met for moderate electromagnetic or gravitational fields. 
Since, according to the equivalence principle, the physical effects of acceleration are locally 
indistinguishable from those of gravity [8], we see that a permanent magnet is unaffected by 
(moderate) acceleration. This conclusion allows us to characterize the accelerating magnet 
by means of equivalent sources in a reference frame S' (with cylindrical space coordinates 
//, 0', z') in which it is instantaneously at rest, then transform this to the laboratory frame 
S (with coordinates p,(j),z). In the laboratory frame, the origin of S' is specified by the 
coordinates [0, 0, z M (£)], where z M (t) is the z-coordinate of the center of mass of the magnet 
and ZM{t)z its velocity. Here and throughout, a caret denotes a unit vector. 

The magnetization vector of an axially symmetric magnet may be represented as 
_M'(p', z') = mP(p', z')z' in its instantaneous rest frame S', where M! is the magnetization, 
mil is the magnetic dipole moment, and P(p', z') is an indefinite density function whose 
integral is normalized to unity. It is important to realize that in general the magnet is in 
accelerated motion, so that this equation embodies the stipulation that the magnetization 
of a permanent magnet is unaffected by acceleration. The effective (or "bound") current 
density corresponding to the above magnetization is given by Jm=V x M. [9]. Therefore, 
we find J'm(p', 0', z') = —m[dP(p', z')/dp']<f>' in the rest frame of the magnet. This quantity 
must now be transformed to the laboratory frame. 

At this juncture we recall that the charge and current density together transform as a 4- 
vector under a Lorenz transformation [10] . Here the charge density of the magnet vanishes in 
its rest frame, and since the current density vector is transverse with respect to the direction 
of relative motion, the charge density in the laboratory frame must vanish as well. As a result 
the transformation equations reduce to the simple result that Jm(p, z, (ft, t) = J'm(p', </> ', z'). 
Note the emergence of the time dependence in J M (p,(fi,z,t), which originates in the fact 
that z' depends on time as well as on z. Indeed for nonrelativistic speeds the coordinates 
transform as in p = p', <fi = </>', and z' = z — ZM{t), and we find 

T I A \ dp[p, Z ~ Z M {t)} ~ m 

J M {p, 0, z) = -m 0, (1) 

which expresses the effective current density of the magnet in the laboratory frame S. We 
note here that the restriction to nonrelativistic speeds in Eq. (1) can be removed by simply 
setting z' = ^[z — z M (£)], where 7=[1 — v 2 /c 2 ] and v = z M (t). 

The case of a uniformly magnetized cylindrical magnet which will be studied in detail 
later corresponds to P(p, z) = 0(~ — \z\)Q(a — p) / (iia 2 L) , where a and L are the radius and 
length of the cylinder respectively and is the standard step function. Using Eq. (1), we 



find 

3tT(p, 0, z, t) = -^Q^ -\z- z M (t)\)5(a - p)4> (2) 

7ra L Z 

for the effective current density of a uniformly magnetized cylindrical magnet in the labo- 
ratory frame. 

B. Vector Potential of the Moving Magnet 

Our task here is the calculation of the vector potential corresponding to the current 
density distribution of the moving magnet given in Eq. (1). Since this current density is 
transverse (or divergenceless) and the charge density vanishes, the Lorenz and Coulomb 
gauges are equivalent here, with the common gauge condition given by V • Am = 0, where 
Am is the vector potential of the moving magnet. The electromagnetic fields of the magnet 
are obtained from B = V x A and E = —dA/dt. The azimuthal symmetry of the current 
density, on the other hand, allows us to write A M (p,(j),z,t) = A M (p,z,t)<f). Thus, like 
the current density, the vector potential has an azimuthal (</>) component only that does 
not depend on the azimuthal coordinate 0. This implies that (a) the electric field is also 
purely azimuthal, and (b) the magnetic field has a radial (p) as well as a longitudinal (z) 
component. It is the radial component of the magnetic field arising from sources external 
to the magnet that exerts the drag force on the moving magnet. 

Using the standard solution for the vector potential in the quasi-static limit [12], we find 

A M (p, z,t)= £ f dz'd<f)'p'dp'[p 2 + p' 2 - 2pp' cos(0 - 0') + {z- z'ff 1 * 

xJ M (p\ z\ £)cos(0-0'), (3) 

where we have used </> • <fi' = cos(0 — (/)'). 

Next we substitute Jm from Eq. (1) into Eq. (3) while making use of the representation 

[13] 

— - 2 f + °° 
[p 2 + p' 2 - 2pp' cos(0 - 0') + {z- z') 2 } ~ 2 = - / dk exp[ik(z - z')\ 

x{ l -I (\k\p')K (\k\p) + ^ cosK0 - <j>')]I n {\k\p')K n {\k\p)}, (4) 

which is valid for p > p' . For < p < p f , p and p' must be interchanged in Eq. (4). 
The above substitution leads to 

A M (p, z,t) = -^ I ' dz'dp'dkex V {ik[z - z M (t) - z'}} 9 ^'^ p'h{\k\p')K x {\k\p), (5) 
Ztt J op 

where the fc-integral ranges over (— oo, +oo) as before. A more convenient form of this 
equation obtains if we use the Fourier representations 

/+oo 
dkexp{ikz}p(p,k) (6) 

— oo 

and 

^•+00 

A M {p,z,t) = {2n)- 1/2 dkexp{ik[z-z M {t)]}A M {p,k). (7) 



Using these representations, we find from Eq. (5) 

A M (p, k) = -/i m f™ dp M !!) k) P 'h{\k\p)Ki{\k\p)- (8) 

J o dp 1 

We recall here that for < p < p 1 , p and p' must be interchanged in Eq. (8). 

It will prove convenient to define a function b (k) by setting A M (p, k) = bo(k)Ki(\k\p). 
Thus we have from Eq. (8), 

b (k) = -pomj^dp'^^p'hilklp'). (9) 

For the case of a uniformly magnetized cylinder, we find from Eqs. (1), (2), (6), and (9), 

6 ° (fc)= w (27r) ~WJW /l(l kl a) - (10) 

III. ELECTROMAGNETIC FIELDS OF THE MAGNET-PIPE SYSTEM 

Having developed the fields of the moving magnet in free space, we now turn to finding 
the fields of the magnet-pipe system. We do this by developing the general solution followed 
by the imposition of continuity conditions. 

A. General Solution 

Our first task here is finding the governing equations for the vector potential in the three 
regions (i) a < p < R X) (ii) R\ < p < R 2 , and (iii) R 2 < p, where R x and R 2 are the inner 
and outer radii of the pipe, respectively [16]. We will actually treat the case of a medium 
with permeability p [17] and conductivity a corresponding to region (ii). For regions (i) and 
(iii), we will simply replace p and a with p Q and 0, respectively. Recall from our discussion 
in section I B that the fields of the magnet-pipe system can be safely calculated in the quasi- 
static limit where displacement currents are neglected. Therefore the equation obeyed by 
the vector potential in the Lorenz (or radiation) gauge reduces to 

(V 2 - pad/dt)A = 0. (11) 

As we saw in section II B, the azimuthal symmetry of the system allows us to represent the 
vector potential in the form A(p, <p, z, t) = A(p,z,t)4>. Moreover, A(p,z,t) is conveniently 
represented as a Fourier integral following the example of Eq. (7): 

/+00 
dkexp{ik[z-z M {t)]}A{p,k). (12) 

—00 

Using this representation in Eq. (11), we find 

(FT + -7j 2 ~ fc2 + rtpiav)A(p, k) = 0. (13) 

dp p op p z 



The general solution of Eq. (13) is a linear combination of I\{y K 2 p) and Ki(yn 2 p), where 
k 2 = k 2 — ipavk [18]. Here and below we will use y/£ to denote that root of £ which has a 
positive real part when £ does. Of course when £ is real and positive, this notation reduces 
to the standard convention whereby \/£ stands for the positive root of £. Note also that 
K\{yK?p) is singular at p = and vanishes exponentially as p — > oo, while Ix{yK*p) vanishes 
at p = but diverges exponentially as p — > oo. Using this information, we construct the 
solution to Eq. (13) in the three regions as follows: 

A«(p, k) = A M (p, k) + h.{k)h{\k\p), (14) 

AM{p, k) = hiQK^Vtfp) + h(k)h(Vtfp), (15) 

A^\p,k) = b 4 (k)K 1 (\k\p), (16) 

where, it may be recalled, A M (p,k) = b (k)Ki(\k\p) is the term corresponding to_the po- 
tential of the moving magnet in free space and bo(k) is given in Eq. (9). Above, A^ n \p, k) 
represents the solution of Eq. (13) in the nth region. Note that we have set a = for regions 
(i) and (iii) as stipulated. 

Using the analogy of waves, one may interpret the terms appearing in Eqs. (14-16) as 
"reflections" and "transmissions" resulting from the "incident" term Am(p, k). This source 
term representing the contribution of the moving magnet is incident upon the inner surface 
of the pipe. The second term in region (i) is the reflection from the inner surface of the 
pipe into the interior. The two terms in region (ii) correspond to a linear combination of 
the transmitted term from region (i) and the reflected term from the outer surface of the 
pipe. In region (iii) we only have the transmitted term from region (ii), since there will 
be no reflection from "the surface at infinity." From a mathematical point of view, on the 
other hand, one starts with a linear combination of the solutions of Eq. (13) in each region 
and proceeds to impose the required boundary conditions. Thus in region (i), the singular 
term [i^i(|fc|p), singular at p = 0] is normalized to represent the contribution of the moving 
magnet, while in region (iii), the singular term [^(j Arjp), singular at p = oo] is excluded 
to ensure that the fields vanish far from the magnet-pipe system. Of course the physical 
sources of all six terms are the (bound) magnetization currents in the magnet and in the 
pipe wall as well as the conduction currents in the pipe wall. 



B. Continuity Conditions 

Our next task is the formulation of continuity conditions across the two boundary surfaces, 
the inner and outer surfaces of the pipe wall. We recall from above that the electric field 
is purely azimuthal, while the magnetic field has radial and longitudenal components. We 
will denote these E^ (= <fi ■ E), B p (= p ■ B), and B z (= z ■ B), respectively. Now the 
boundary conditions require the continuity of Eq (Faraday's law), B p (absence of magnetic 
monopoles), and H z (the Ampere-Maxwell law and absence of surface currents) across the 
two boundary surfaces, where H = B/u in regions (i) and (iii), and H = B/u in region 
(ii). When expressed in terms of A, the first two of these conditions require the continuity 
of A while the third condition requires the continuity of fi~ 1 d(pA)/dp. An inspection of 
Eq. (10) shows that these conditions must also be obeyed by A. This gives us the continuity 



equations that must be imposed on the solutions of Eq. (11). Recall that we must also apply 
the conditions k — > k and p — > po in regions (i) and (iii). 

Upon imposing the above-stated continuity conditions on the solutions given in Eqs. (12- 
16) at the boundary surfaces p — R\ and p = R 2 , we find the following set of equations: 



b {k)K 1 {\k\R l ) + b l {k)h{\k\R 1 ) = b 2 (k)K 1 (VK 2 Ri)+ & 3 (fc)Ii(V« 2 iii; 
b 2 {k)K 1 {y/^R 2 ) + b 3 {k)I l {^R 2 ) = b 4 (k)K\(\k\R 2 ), 



^[-boWKoQklRj + &i(fc)io(|A;|i*i)] = —[-b 2 (k)K ( \^R 1 ) + b 3 (k)I (V^R l )}, 

Po p 

2 I 1 I 

b 2 (k)K (^R 2 )] + b 3 (k)I (^R 2 )} = l -^[-b,(k)K (\k\R 2 )}. (17) 



P Po 

This linear set can be solved by standard methods to find the unknown coefficients bi(k) 
through b&(k). To avoid unnecessary writing, we will only record the solution for bi(k), since 
this is the coefficient that will be needed for the calculation of the drag force in section V: 

&i(fc) = {[K (|A;|i? 1 )K (|A;| J R 2 )T 11 +/5Ko(|fc| J Ri)K 1 (|A;|i? 2 )T 1 o 

+ [I Q {\k\R x )K Q {\h\R 2 )T xl + f3I (\k\Ri)Ki(\k\R2)T w 
+pi 1 (\k\R 1 )K Q (\k\R 2 )T Q1 + /5 2 / 1 (|A;|/2 1 )K 1 (|A;| J R 2 )Too]}6o(A;), (18) 



where 



T 00 = K (a\k\Ri)I (a\k\R 2 ) — I (a\k\Ri)K (a\k\R 2 ) 

T i = i^ (tt|fc| J Ri)/i(tt|fc| J R 2 ) + T («|A;| J Ri)i^i(tt|fc| J R 2 ) 

T10 = K^alklR^IoialklR-z) + hialklR^Ko^alklR^ 

T n = K^alklRjI^alkllh) - Ii(a\k\Ri)Ki(a\k\Rz), 



(19) 



and 



a - Vny\k\ l l _ l ^!^l. 



P=^ = ^-Jl- l f^^. (20) 

prel Prel V k 

Here p re \ = p/po represents the relative permeability of the pipe material. It is worth noting 
that a and f3 have positive real parts by construction. 

With the coefficients given in Eqs. (18-20), we have in Eqs. (9), (12), and (14-16) a 
complete solution to the magnet-pipe system in all regions. 



IV. DRAG FORCE ON THE MOVING MAGNET 

We are now in position to calculate the braking force exerted on the moving magnet. 
Recall that the density of magnetic force exerted at a point r of a current distribution J(r) 



is given by J(r) x B(r), where B(r) is the magnetic field strength at the point in question. 
Thus the force exerted on the moving magnet is given by 

d 3 r J M (r,t)xB(r), (21) 

where Jm( t , t) is the effective current density of the magnet given in Eq. (1). Note that B(r) 
can be replaced with B exi (r) in Eq. (21), where the latter is the magnetic field produced in 
region (i) by sources external to the magnet. Using Eqs. (12) and (14), we find 

/+oo 
dkexp{ik[z - z M {t)]}bi{k)I x {\k\ P )j>. (22) 

— oo 

Combining Eqs. (1), (6), and (22), and performing a few straightforward operations, we find 
from Eq. (21) 

F = -2mmz[ kdk [ pdp6i(fc)/i(|fc|p) ^^' ~^ . (23) 

J -oo JO Op 

At this juncture we use Eq. (9) to rewrite this result in the form 

/+oo 
kdkb (-k)bx(k). (24) 

— oo 

The quantity b\(k) was defined in Eqs. (18-20) et seq, and can conveniently be expressed 
as b\{k) = Q(k)bo(k). An inspection of Eqs. (18-20) shows that the real and imaginary parts 
of Q(k) so defined are even and odd functions of k, respectively. On the other hand, Eqs. (9) 
and (6) show that bo(k) = b^—k). These properties allow us to rewrite Eq. (24) as 

/+oo 
fcdfc|&o(fc)| 2 Im[Q(fc)], (25) 

o 

where we have also replaced z by v to make the retarding nature of the force manifest. This 
equation gives the drag force acting on the moving magnet and is a central result of our 
analysis. 

Using Eq. (10) which gives 6g m (fc) corresponding to the uniformly magnetized magnet, 
we find from Eq. (25) 



sin(A;L/2) 
2tt 2 J |_ (kL/2) 



v^- / dkk 



Ii(ka) 
Jka/2) 



Im[Q(fc)]. (26) 



This equation gives the drag force acting on the moving magnet for the case of a uniformly 
magnetized cylinder. 

It is instructive to rederive the drag force formula in Eq. (25) from energy conservation. 
Let us consider the magnet-pipe system at a moment the magnet is moving through the pipe 
with velocity v. Under the quasi-static conditions stipulated in the section IB, ohmic power 
dissipation in the pipe wall must be balanced by a decrease in the kinetic energy of the moving 
magnet, since the power flow into the electromagnetic field configuration, including radiation, 
is negligible. But the decrease in kinetic energy corresponds to a resistive force according 
to the work-energy theorem, P h m = — F ■ v, where P hm is the rate of ohmic dissipation (or 

9 



Joule heating). Needless to say, the magnet may be experiencing non-electromagnetic forces 
such as gravity or air drag which would have their own power contributions. Now Poynting's 
theorem assures us that the dissipated power equals the flux of Poynting's vector into the 
pipe wall through its inner surface. Putting these two observations together, we arrive at 



F • v = - / pd(j)dzS ■ p, (27) 

where S is Poynting's vector, given by = /i ( ^ 1 E x B. Since by symmetry F must have an 
axial direction, we can rewrite Eq. (27) as 



, 1 f 

-v / pd(j)dzp 

VovJ p=Rl 



§ x (V x A) 



(28) 



where we have used B = V x A and E = —dA/dt. 

At this juncture we use Eq. (12) to substitute the Fourier representation of A- % > into 
Eq. (28). This allows us to perform all implied integrations except one, with the result 



/oo 
kdk 
— oo 



A«(p,k) eAW £-» 



(29) 

P=Ri 



Next, we use Eq. (14) to replace A^\p,k) with its expression in terms of modified Bessel 
functions. Of the resulting four terms inside the square brackets, two are even in k and make 
no contribution to the integral in Eq. (29). The other two terms equal \k\b Q (k)bi(— k)W, 
where W = l/\k\R\ is the Wronskian of Ki(\k\Ri) and Ii(\k\Ri) [18]. Upon replacing bi(—k) 
with b (—k)Q(—k) in Eq. (29), and recalling that the real and imaginary parts of Q(k) are 
even and odd functions of k respectively, we recover Eq. (25). This completes the derivation 
of the magnetic drag force from energy conservation. 

V. PROPERTIES OF THE DRAG FORCE 

The remainder of this paper is devoted to a detailed discussion of the properties of the 
drag force for the case of a uniformly magnetized cylinder, given in Eq. (26). It should be 
remembered, however, that the magnetization distribution of a typical magnet depends on 
its type and manufacturing method, and almost certainly deviates from uniformity. The 
following study is thus intended to provide a benchmark that typifies general properties. 

We shall begin our study by characterizing the main features of the drag force here, 
including its dependence on the shape and speed of the magnet as well as the material 
properties of the pipe wall. In the appendices, we will derive asymptotic expressions for the 
drag force of Eq. (26) in a number of physically interesting limiting cases. 

A. Dependence on Magnet Shape 

Let us start by examining the dependence of the drag force on the dimensions of the 
magnet. We have already arranged Eq. (26) in such a way as to isolate and highlight the 
dependence of the drag force on the relevant parameters. The force is opposite the velocity 

10 



[19] and it is scaled by the square of the magnet's dipole moment [20]. For a fixed value of 
the dipole moment, the dependence of the force on the dimensions of the magnet is entirely 
contained in the two bracketed factors in Eq. (26), the first of which depends on the magnet 
length L and the second on its radius a. Each of these factors has been normalized to unity 
at the point-dipole limit, the first at L = and the second at a = 0. An inspection of the 
second factor shows that it increases monotonically with increasing a, which is expected as 
such an increase brings the source currents in the magnet closer to the eddy currents in the 
pipe wall thus increasing the interaction force. On the other hand, the envelope of the first 
factor clearly decreases with increasing magnet length, suggesting a corresponding decrease 
in the drag force with increasing magnet length. This is in fact confirmed by our numerical 
results, and reflects the weakening of the external magnetic field with increasing magnet 
length (with a fixed magnetic dipole moment as stipulated), which field eventually vanishes 
as L/a — ► oo, just as it would for an ideal solenoid. It must be remembered, however, that a 
can at most equal Ri, and L must remain small compared to the distance from the magnet 
to the pipe ends. It is convenient in this regard to define an ordered pair of dimensionless 
parameters, (L/2a,a/Ri), characterizing the inverse aspect ratio of the magnet and the 
tightness of its fit inside the pipe, respectively. Thus for a given dipole strength, the pair 
(0, 1.00) corresponds to a wafer-shaped magnet that would just fit inside the pipe while 
(1, 0.60) describes a "square cylinder" filling 36% of the interior cross section of the pipe. 

Numerical results were obtained from Eq. (26) using a Mathematica program developed 
for this purpose [7]. Figure 1 shows a plot of the drag force versus magnet speed for m = 1.00 
A m 2 , fi rel = 1.00, a = 5.00 x 10 7 ft" 1 m -1 , R x = 10.0 mm, R 2 = 11.0 mm, and five different 
shapes, (a) (typical cylinder, loose fit) =>- (|,0.60), (b) ("square" cylinder, loose fit) =>- 
(j,0.60), (c) point-like cylinder =>- (y,^ 0), (d) (short cylinder, snug fit) =>- (|, 0.96) (e) 
(circular wafer, loose fit) =^ (~ 0,0.60). The dashed line, on the other hand, represents the 
idealized model limit derived in appendix B, F tdl = — v45/i 2 1 m 2 w<7,s/1024R 1 4 , with the pipe 
thickness s set equal to R 2 — Ri = 1.0 mm in order to facilitate comparison to the exact 
results. 

The five cases shown in Fig. 1 demonstrate the effect of the shape of the magnet on 
the drag force. A comparison of cases (a) and (b) in Fig. 1 clearly shows that for fixed 
values of the dipole moment and speed, the drag force increases as the magnet is shortened. 
Comparing cases (b) and (d), or (c) and (e), on the other hand, we see an increase in 
the drag force with increasing magnet diameter. These five cases clearly demonstrate the 
strong influence of the shape of the magnet on the drag force. In particular, they clearly 
demonstrate the quantitative inadequacy of the point-dipole approximation, case (c), across 
a broad range of speeds, as discussed in section I A. 

The dashed line representing the idealized model corresponds to augmenting the point- 
dipole approximation by the thin- wall assumption, (R 2 — R±)/Ri — > [21]. A comparison of 
the dashed line and curve (c) in Fig. 1 shows a relative deviation of about 10% in the low- 
speed regime. This deviation should be compared with the ratio (R 2 — Ri)/R\ characterizing 
the relative thickness of the pipe, which equals 0.100 for all cases in Fig. 1. It is clear that the 
error caused by the thin-wall assumption is, in general, unacceptably large for a reasonable 
agreement with measurement results under typical conditions. Exceptions can of course 
occur: for case (d) corresponding to a short cylindrical magnet fitting snugly inside the 
pipe, the errors caused by the point-dipole and thin-wall approximations nearly cancel one 
another in the low-speed regime. 

It is appropriate at this juncture to compare the prediction of our analysis to the measured 
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FIG. 1: Plot of the drag force versus the magnet speed for a fixed value of the dipole moment and 
four different shape parameters (L/2a,a/R\) (see the text for definitions): (a) typical cylinder, 
(|,0.60), (b) "square" cylinder, (y,0.60), (c) "point-like cylinder" (y,— 0), (d) short cylinder 



(|, 0.96), and (e) circular wafer, (~ 0, |). The dashed line represents the idealized limit. 



results of MacLatchy et al. [3], bearing in mind the important caveat that the magnetization 
distribution in their experiment was not uniform. Using the parameter values reported by 
these authors, namely m = 0.67 A m 2 , a = 5.08 x 10 7 Q~ x m" 1 , R± = 7.29 mm, R 2 = 7.96 
mm, L = 6.4 mm, and a = 6.3 mm, and equating their reported magnet (plus tape) weight 
of 0.060 N to the magnetic drag force of Eq. (26) as well as setting /i re / = 1.00 appropriate 
for copper, we find a magnet speed of 11.9 ± 0.5 cm s _1 . MacLatchy et al. reported a 
measured terminal speed of 12.7 ± 0.4 cm s _1 , and compared this to the prediction of the 
idealized model, 17.8 cm s _1 . Note that in this comparison we are neglecting air drag on 
the magnet, estimated by the authors to be less than 0.1% of the magnetic drag force hence 
deemed negligible. The quoted error of 0.5 cm s _1 in our calculated result is estimated on 
the basis of the (implied) precision level of the measured values given by MacLatchy et al. 
for L, a, m, R 2 — -Ri, and the magnet weight. Keeping in mind that the magnetization 
distribution of the "button" magnet used in the experiment was not uniform, and the fact 
that the drag force is rather sensitive to this distribution as demonstrated by the numerical 
modeling of MacLatchy et al, we conclude that our predicted value of the magnet speed 
agrees with the measured value within the uncertainties. 

The circular wafer (L = 0) geometry merits special attention since, in that case, the drag 
force grows without limit as a —* R\. In other words, the shape characterized by (0, 1) is 
a singular limit where the magnetic drag force diverges. This behavior may be understood 
as follows. In the limit of (0, 1) geometry, the effective (or "bound") source currents of the 
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magnet are concentrated in a single current ring of radius R±. However, such a source would 
induce similarly concentrated eddy current rings in the pipe wall on its interior surface, 
i.e., at zero distance from the source itself, thus causing infinite interaction forces due to 
overlapping current rings. This situation is analogous to the divergence of the image force 
on a point charge as it approaches the surface of a conductor. It is worth remarking that 
while the limit itself is unphysical, the sharp increase in the magnetic drag force for L <C a 
and a — > R\ is very real and already discernible from a comparison of cases (c) and (e) in 
Fig. 1. In practice, this sharp increase in the magnetic drag would be accompanied by a 
parallel increase in the air drag force caused by the unavoidable constriction in the flow of 
air around the magnet [22]. 

The effect just discussed can also be understood in reference to Eq. (26). Using the asymp- 
totic properties of the modified Bessel functions [18], we find that Im[Q(A;)] — > exp(—2kRi)/k 
as k — » oo. Now for L/a — > and a —* R\, the shape factors in Eq. (26) behave like 
exp(2kR 1 )/k 3 as k — > oo. Putting these two statements together, we conclude that the inte- 
grand in Eq. (26) behaves like k~ l for large k, which implies a logarithmic divergence at the 
upper limit of the integral. Recalling that k is the Fourier conjugate of z [cf. Eq. (26)], so 
that large values of k are associated with short values of z, we conclude that this divergence 
is a short-distance effect. This is of course the conclusion we reached above on physical 
grounds. 



B. Dependence on Magnet Speed 

The dependence of the drag force on the speed of the magnet is contained in Q(k) 
through the combination fiav and is more involved than the dependence on its shape. It 
proves expedient to discuss this dependence in terms of a length parameter defined by 
^o = (fJ'O^v) - . As discussed below, Iq represents the penetration depth of the fields into the 
pipe wall under appropriate conditions. The two distinguished ranges of v, which we will 
refer to as "low" and "high" speed [23], correspond to l ^> R\ and Iq <C R±, respectively 
[24] . We will first consider the low-speed regime, as it is the relevant one in practice. In this 
speed range, we have a ~ 1 — i/(fi~^l k), so that akR 1 ~ kRi — i(Ri//j,~^l ), and similarly 
for akR^. Thus the imaginary part of the argument of various modified Bessel functions 
is much smaller than the corresponding real part. This would in turn imply that, to the 
leading order, the imaginary part of those functions is linear in Ri/Iq and i^Ao- This fact 
implies the same for Q, to wit, that in the low-speed regime the imaginary part of Q is linear 
in Ri/Iq (and R2/I0), hence proportional to o~v [25]. This confirms the expectation that in 
the low-speed regime where Iq 3> R±, the braking force is a linear drag proportional to the 
conductivity of the pipe wall. 

The linear nature of the drag force in the low-speed regime is in evidence for all cases in 
Fig. 1. In typical demonstration setups, one would expect terminal speeds of the order of 
1 ms" 1 or less, which corresponds to the neighborhood of the origin in Fig. 1 where the drag 
force is closely proportional to the magnet speed. Indeed from the the low-speed condition 
ficrvRi <C 1 given above, we can estimate the relative deviation from linearity to be of the 
order of (/xavRi) . For case (d) of Fig. 1, this estimate of deviation gives (0.63?; [ms -1 ]) , 
which amounts to about 2.5% at a terminal speed of 25 m s _1 . 

In the high-speed regime (which would require the projection of the magnet into the pipe 

with a suitably high initial speed) where Iq <C Ri, we have a ~ (1 — ?)/(2/i~ e 1 /o^) 2 for values 
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of kR\ ~ 1 which provide the main contribution to the integral in Eq. (26). This implies 

that akRi ~ (1 — i)(kRi/2)^ (R 1 /fx~^ l l ) 2 , and similarly for akR 2 . Since Ri/l ,R 2 /l 3> 1 
in this speed range, the arguments of the modified Bessel functions in Q which involve a 

will scale as (Ri/Iq) 2 or (R 2 /lo) 2 f° r ^1 ~ 1- Since these functions behave exponentially 
for large values of their argument, the said scaling behavior results in an overall suppression 
Im(Q) in the high speed regime, corresponding to a decrease in the drag force. In other 
words, the drag force is a decreasing function of av in the high-speed regime. Physically, a 
pronounced skin effect which suppresses the penetration of the field into the pipe wall takes 
hold at high speeds, thereby reducing the eddy currents and the drag force resulting from 
them. Indeed, recalling from section IB that the dominant time scale in the magnet-pipe 
system is Ri/v, which corresponds to a frequency u> ~ v/Ri, we see that the high speed 

condition l <C R 1 is equivalent to the inequality (/i<ru;)~ 5 <C R\ [24]. But this last condition 

states that the skin depth corresponding to the dominant time scale is much smaller than 

Ri and, barring unusually thin-walled pipes, much smaller than R 2 — R\ as well. In other 

_ 1 
words, we have (fiau) 2 <C R 2 — R\ in this limit, which is the condition for a pronounced 

skin effect [5] . It is appropriate to repeat here that the high speed limit is not likely to occur 

in typical arrangements of the magnet-pipe demonstration. 

The nonlinear behavior of the drag force at intermediate speeds and its eventual drop 

at high speeds deduced above are clearly in evidence for all cases displayed in Fig. 1. Note 

that the high-speed decline in the drag force becomes sharper as the source currents in the 

magnet are more highly concentrated. Recalling our discussion of the singular geometry 

above, we see the reason for this behavior: as source currents become more concentrated, 

short- distance, equivalently high-Zc, contributions become more important, a feature that 

runs contrary to the high-speed situation where high k values are relatively less important. 

C. Dependence on Conductivity and Susceptibility 

Recall that the drag force dependence on the magnet speed is through the combination 
fxav. This implies that the drag force behavior versus a follows the same pattern as for 
v. In particular, higher conductivity makes for a larger drag force at low speeds hence the 
use of copper tubes for demonstration purposes. At high speeds, however, the drag force 
decreases with pipe conductivity, just as it does with magnet speed. An interesting conclu- 
sion, therefore, is that the drag force vanishes for a perfect conductor which, as explained 
above, is simply a consequence of a strong skin effect. This behavior is displayed in Fig. 2, 
which is a plot of the drag force versus a for case (d) of section VA at a magnet speed 
of 0.10 m s _1 . It is important to note that the decline of the drag force for high values 
of a is not directly applicable to practical arrangements which usually correspond to the 
low-speed (or possibly intermediate-speed) regime . Also to be noted is the fact that the 
designation "perfect conductor" here implies electric conduction without resistance, and is 
to be distinguished from "superconductor" which in addition implies a distinct magnetic 
behavior as noted below. Of course inasmuch as a superconductor has zero resistance, the 
above argument shows that the magnetic drag force vanishes for a superconducting pipe. 
Needless to say, this conclusion as well as the one above for a perfect conductor, directly 
follow from the energy conservation principle, which in this case asserts that there can be 
no drag force without a corresponding dissipation of power in the pipe. 

This brings us to the question of how the magnetic properties of the pipe wall affect the 
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FIG. 2: Plot of the drag force versus the conductivity of the pipe for case (d) with v = 0.10 ms 1 



drag force [17]. Thus far we have assumed a linear magnetic material of relative permeability 
fi re i for the pipe wall. Inasmuch as the magnetic susceptibility of non-ferromagnetic common 
metals differs little from that of the vacuum, one can set \x re i ~ 1 for practical purposes, as we 
did for the cases displayed in Fig. 1. However, it is physically interesting and meaningful to 
consider the extreme cases of /i re ; — > and fi re i — > oo, corresponding to perfect diamagnetism 
and perfect paramagnetism, respectively. Perfect diamagnetic behavior is exemplified by a 
(type I) superconductor which would exclude any magnetic field from its interior (save for 
a very small penetration depth). This phenomenon, known as the Meissner effect, is not a 
mere consequence of perfect conductivity and serves to distinguish a superconductor from a 
material that conducts electricity without dissipation [26]. Perfect paramagnetism, on the 
other hand, is approximately realized by "soft" ferromagnetic materials which have a narrow 
hysteresis loop and can be idealized as highly susceptible, linear magnetic materials. 

An inspection of Eq. (26) shows that relative permeability occurs in the quantities a 
and P = a/fi r ei through the combination fi re ifi av in a, and in addition, j3 is inversely 
proportional to [i re \. Leaving the latter aside for the moment, we conclude that the drag 
force dependence on /j, rel is much like its dependence on the magnet speed. Thus the drag 
force must vanish for both fi re i —>■ and \i re x —>■ oo, corresponding to strongly diamagnetic 
and paramagnetic limits. As further discussed in appendix D, the additional dependence 
of j3 on fi re i does not change these qualitative features. Thus the drag force corresponding 
to a superconductor vanishes not only because of its perfect conductivity but also due 
to its perfect diamagnetism. Similarly, the drag force is seen to be small for a strongly 
paramagnetic conductor, or soft ferromagnetic alloys that behave like one. The asymptotic 
behavior of the drag force for \i re i <C 1 and \i re \ ^> 1 is analyzed in appendices D and E, 
respectively, where the foregoing conclusions are explicitly confirmed. 

Figure 3 shows a plot of the drag force versus /i re / for case (d) of section V A at a magnet 
speed of 1.0 m s _1 , where the features deduced above are in evidence. 

VI. CONCLUDING REMARKS 

The treatment of the magnet-pipe system in this paper has been based on axial symme- 
try. Therefore, any deviation of the magnet from a linear, axially centered motion such as 
wobbling or tumbling would violate this underlying assumption and cause a departure from 
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FIG. 3: Plot of the drag force versus the relative permeability of the pipe for case (d) with v = 1.0 
-l 
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the predicted results [25] . Where necessary, the magnet can be placed inside an electromag- 
netically inactive casing with an optimum profile for stability. This procedure would also 
serve to maintain a fixed air drag coefficient when comparing magnets of different profile. 
Similarly, the analysis in this paper assumes an infinitely long pipe, so that in practice the 
magnet ends must be many times the interior pipe diameter away from the pipe ends when 
measurements are taken. 

The results of this paper can be used for precision magnetic braking studies and appli- 
cations. The computer program provided for use with this paper [7] has been written for 
the case of uniform magnetization, but can readily be modified to deal with the general case 
using Eq. (25). For precision studies, it may be preferable to use an electromagnet instead of 
a permanent magnet, since the former allows a more convenient characterization of sources 
and fields. 
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APPENDIX A: LOW MAGNET SPEED 

Here and in the following we will outline the development of a few limiting expressions 
for the drag force given in Eq. (26). The methods used are those of approximation and 
asymptotic analysis of integrals, the details of which would take us beyond the scope of this 
paper [27]. 

As a preliminary step we recall that in all cases except for appendix D the main contri- 
butions to the integral of Eq. (26) originate from the region k ~ R^ . This fact makes it 
convenient to rescale the integration variable therein according to k — u/R\. Effecting this 
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substitution, we obtain 

^ r+x r sm(uL/2Rx 



'-^%l dun* 



(uL/2Rx 



I\[uajR\ 
_ {ua/2R 1 



2 

lm[Q(u, s/Rx, Ri/lo, fi re i)}, 



(Al) 

where we have explicitly defined the dependence of Q on three dimensionless parameters 
which characterize the dimensions and material properties of the pipe as well as the magnet 
speed. Note that s stands for the thickness of the pipe wall here. We shall use this rep- 
resentation to derive asymptotic expressions for the drag force in a number of interesting 
limiting cases. 

The low-speed regime is defined by the condition Rx/Iq <C 1 [24]. To establish the fact 
that F um is linear in v in the low-speed limit, as discussed in Sec. VB, we note that 
Q(u, s/Ri, 0, fi re i) is real at v = (corresponding to the vanishing of the drag force at zero 
speed), so that an expansion of Im(Q) about v = using the asymptotic properties of the 
modified Bessel functions [18] yields Im[Q(u, s/Rx, Ri/lo, Hrei)] — £Rx/Iq f° r small Rx/Iq- 
The quantity C equals lm[dQ(u, s/Rx, 0, jJ, re i)/d(Rx/lo)], and is given by a long expression 
which need not be reproduced here. Substituting the approximate form of Im(Q) in Eq. (Al), 
we immediately obtain the linear drag behavior in the low-speed limit: 

F is P ^ _ Cavir [n rel ^ av R x «!), (A2) 

where C depends on all parameters other than a and v. 

APPENDIX B: THE IDEALIZED MODEL 

The idealized model involves three assumptions: (a) the point-dipole limit, L/Rx — > 
0, a/Rx — > 0, (b) the low-speed approximation, Rx/Iq <^ 1, and (c) the thin- wall approxi- 
mation, s/Rx ^ 1. The point-dipole limit is readily implemented by setting the two shape 
factors equal to unity. Items (b) and (c), on the other hand, require applying the thin- 
wall approximation to the quantity C introduced in Eq. (A2). Using the properties of the 
modified Bessel functions, we find from Eqs. (18-20) that C — > (s / Rx)u[Kx(u)} 2 . Therefore, 

Q idl = iii Q {s/Rx)u[Kx{u)fav, (Bl) 

where Q ldl represents Q in the ideal limit. Substituting this result in Eq. (Al) (with the 
shape factors set to unity), we find 

— - ^"^v, (B2) 



1024/V 



in agreement with previous results [2]. 



APPENDIX C: HIGH MAGNET SPEED 

As stated in Sec. VB, in the high-speed regime where Iq/Rx "C 1, the quantity a in 
Eqs. (18-20) has a large magnitude which forces the respective modified Bessel functions to 
their asymptotic range. Since the four quantities Ty in Eqs. (19) have the same asymptotic 
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limit [18], they cancel out of Eq. (18). The remaining terms can then be simplified using 
the properties of the modified Bessel functions, leaving the result 

Im (Q) = i^ .^inn^Lit-ip^ GWo<™ » Ri). (CI) 

kRi [I (\k\Ri)] + l^lhdklRj] 

The remaining v dependence in Eq. CI resides in /3 — a/fi re i. Recalling from Sec. VB that 
a — > (1 — i)/(2fj,~ el lok) 2 in the high-speed limit, we can reduce the above equation to [24] 



Im(Q) ^ v^ 

V2fc// ^i2i[/i(|fc| J Ri)] 2 



—j (HrelHoO-V > i?i). (C2) 



^uEi/Zof/i^)] 
Upon substituting this expression in Eq. Al, we find 

F hsp = -^^^(a/i^L/i^/^v, (C3) 

where T\ is a form factor which depends on the dimensions of the magnet as fractions of 
the interior diameter of the pipe. It is defined by 

^(a/R 1 ,L/R 1 )= g ^^\ (C4) 



where 



/+oo 
duu 5/2 




sin(nL/2i?i 

(uL/2R 1 ) 



I\(ua/R\) 
_ {ua/2R 1 ) 



[h{u)Y 2 . (C5) 



Note that the form factor T\ has been normalized to unity at the point-dipole limit. 

The asymptotic formula in Eq. (C3) is valid in the high-speed limit / u re //i crt' ^> R\ and 
describes the behavior of the drag force for very high magnet speed or pipe wall conductivity 

[24]. 

APPENDIX D: HIGHLY DIAMAGNETIC PIPE 

This is the limit fi re i — * 0, which is appropriate to a highly diamagnetic pipe wall ma- 
terial. In this limit the factor /3, which occurs in Q, grows large and severely suppresses 
the magnitude of Im(Q) in Eq. (Al), and in addition limits its contributions to very small 
values of u. This mathematical behavior reflects the physical phenomenon of magnetic flux 
expulsion that accompanies the approach to perfect diamagnetism. 

Since only small values of u are important in Eq. (Al), we can replace all modified Bessel 
functions by their asymptotic values in the small argument limit [18]. This leads to the 
following approximate expression for Q: 



a 2 \n(R 2 /R 1 
fi re i + lu 2 ln(R 2 /Ri)a 2 



q = - .r™,. -(-<i) Pi) 
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Next, we find Im(Q) and change the variable of integration in Eq. (Al) according to w = 

-i i 

// r Jln(i? 2 /-Ri) 2 u. After some calculation and simplification, we find 

/+oo i —2 

ciCC 2 [l + 2 C2] *> (D2) 

where 

K = ^^M**/RiTirf2- (D3) 

Finally, we carry out the integral in Eq. (D2) to arrive at the asymptotic behavior of the 
drag force: 

F Mm = -^^[ln(i? 2 /^)]-^ 2 v. (D4) 

The above formula is valid in the limit fi re i — » 0, and is appropriate to a highly diamag- 
netic pipe. 



APPENDIX E: HIGHLY PARAMAGNETIC PIPE 

Here we are considering the limit fi re i — > oo appropriate to a highly susceptible pipe 

wall material [17]. Recall from VC that the dependence of Q on fi re i occurs through the 

combination fi re ifio<jv in a, and in addition through (3 = a/fi re i. Consequently, Eq. CI 

which is appropriate in the high-speed regime holds here as well. Recalling further that 

i 
(3 — > (1 — i)/(2fi re il k) 2 in this limit, we find 

MQ) = -WWj=TT ( .2 — (»rel » i). (El) 

At this point we follow the steps subsequent to Eq. (CI) to find 

F*» = - °- 053 ^ m V o(«/Jt„ L/ Rl ),/™v, (E2) 

R{ V » 

where T§ is a form factor which depends on the dimensions of the magnet as fractions of 
the interior diameter of the pipe. It is defined by 

J r o{a/Ri,L/R l ) = , (E3) 

yo(U,(JJ 



where 



/+oo 
duu z ' 2 




sin(uL/2R 1 

(uL/2R 1 ) 



I\(ua/R\) 
_ {ua/2R 1 ) 



[I (u)r 2 . (E4) 



Note that the form factor JF has been normalized to unity at the point-dipole limit. 

The asymptotic formula given in Eq. (El) is valid in the limit /i re / 3> 1 and describes the 
behavior of the drag force for a highly susceptible pipe. 
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